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Assuming an intrinsic 'Band' shape spectrum and an intrinsic energy-independent emission profile we have investigated 
the connection between the evolution of the rest-frame spectral parameters and the spectral lags measured in gamma-ray 
burst (GRB) pulses by using a pulse model. We first focus our attention on the evolution of the peak energy, Eo,p, and 
neglect the effect of the curvature effect. It is found that the evolution of Eo,p alone can produce the observed lags. When 
Eo,p varies from hard to soft only the positive lags can be observed. The negative lags would occur in the case of Eo,p 
varying from soft to hard. When the evolution of Eo^p and the low-energy spectral index ao varying from soft to hard then 
to soft we can find the aforesaid two sorts of lags. We then examine the combined case of the spectral evolution and the 
curvature effect of fireball and find the observed spectral lags would increase. A sample including 15 single pulses whose 
spectral evolution follows hard to soft has been investigated. All the lags of these pulses are positive, which is in good 
agreement with our theoretical predictions. Our analysis shows that only the intrinsic spectral evolution can produce the 
spectral lags and the observed lags should be contributed by the intrinsic spectral evolution and the curvature effect. But it 
is still unclear what cause the spectral evolution. 
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1 Introduction 

' The phenomenon of the observed spectral lag of Gamma- 
, ray burst (GRB) is very common. Since Cheng et al. (1995) 
' first analyzed the spectral lags of GRBs based on BATSE 
] Large Area Detector channel 1 (25-50 keV) and channel 
' 3 (100-300 keV) light curves, several authors have carried 
' out more analysis work on the GRB lags (e.g., Norris et 
■ al., 2000; Wu & Fenimore 2000; Chen et al. 2005; Yi et al. 
; 2006; Peng et al. 2007). 

Several attempts to explain the origin of the spectral lag 
have been put forward. The dominant behavior observed in 
GRBs, positive lags, is that the soft radiation lagging the 
hard, can be produced by several reasons. An intrinsic cool- 
ing of the radiating electrons cause the radiation to domi- 
nate at lower energies. A Compton reflection of a medium 
at a sufficient distance from the initial hard source will also 
cause a positive lag. While the photons emit at relativis- 
tic speeds from different latitudes of a convex surface (the 
curvature effect) will cause the radiation to be delayed and 
softened, thus producing a observed positive lag (e.g., Qin 
2002; Ryde & Petrosian 2002; Qin et al. 2004; Shen et al. 
2005). The negative lag is defined as the hard radiation lag- 
ging behind the soft. However, few attempts are made to ex- 
plain the origin of the negative lag. The possible scenario is 
a hot medium, say a lepton cloud surrounding a cooler emit- 
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ter, which will cause the soft radiation to be upscattered by 
inverse Comptonization. The photons get harder the more 
scatterings they suffer and thus they are more delayed. A 
definitive explanation of the origin of them has not yet been 
given. The origin of the spectral lags are still unclear. 

Kocevski & Liang (2003) pointed out that the observed 
lag is the direct result of spectral evolution. They measured 
the rate of Ep decay ((f>o) for a sample of clean single- 
peaked bursts with measured lag and used this data to pro- 
vide an empirical relation that expresses the GRB lag as 
a function of the burst's spectral evolution rate. While the 
studies of GRB spectra have shown that GRB spectral evo- 
lution is a common phenomenon (e.g.,, Kargatis et al. 1994; 
crider et al. 1997; Band 1997; Kaneko et al. 2006; Butler & 
Kocevski 2007). 

Recently, Shen et al. (2005) tentatively studied the con- 
tribution of curvature effect of fireball to the lag, and the 
resulting lags are very closed to the observed one. Lu et al. 
(2006) also studied the spectral lags more detailed based on 
Doppler effect of fireball (or in some papers, the curvature 
effect) (see, Qin 2002; Qin et al. 2004; Peng et al. 2006 
for more detailed description) and confirmed that the cur- 
vature effect can produce the observed spectral lags. They 
performed more precise calculation with both formulae pre- 
sented by Shen et al. (2005) and Qin et al. (2004). Some 
conclusions obtained by Shen et al. (2005) were ascertained, 
and more complete conclusions on spectral lags resulting 
from curvature effect were obtained. Lu et al. (2006) argued 
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that, as long as the whole fireball surface is concerned, both 
formulas were identical in the case of ultra-relativistic mo- 
tions. Other cases did not study by Shen et al. (2005) were 
investigated by Lu et al. (2006) and some new conclusions 
were drawn. These seem show that the curvature effect can 
indeed produce the spectral lags. In addition, several studies 
show that the curvature effect also plays an important role 
in the early X-ray afterglow of GRBs where the so-called 
softening phenomenon is observed (Qin 2008a, 2008b; Qin 
2009). 

As pointed out above that the fundamental origin of the 
observed lag is the evolution of the GRB spectra and the 
curvature effect can produce the observed lags. Moreover it 
is mentioned clearly that whilst the curvature effect can ac- 
count for the main part of the hardness ratio curves of the 
GRBs concerned and an intrinsic spectral evolution is re- 
quired to fully explain the observed data (Qin et al. 2006). 
These motivate our investigations below. How the influences 
of intrinsic spectral evolution on the spectral lags are when 
the curvature effect plays an important role or does not? Can 
the negative lags be produced in the course of intrinsic spec- 
tral evolution? We introduce the theoretical formula of the 
GRB pulse in Section 2. Various possible cases of produc- 
ing the spectral lags when neglect the curvature effect are 
studied in Section 3. In Section 4 we investigate the lags 
when spectral evolution and curvature effect take effect si- 
multaneously. We employ a sample to test the prediction in 
Section 5. In the last section, we give the conclusions and 
discussion. 

2 Theoretical formula 

The observed gamma-ray pulses are believed to be produced 
in a relativistically expanding and coUimated fireball be- 
cause of the large energies and the short time-scales in- 
volved. To account for the observed pulses and spectra of 
bursts, the Doppler effect (or curvature effect in some pa- 
per) over the whole fireball surface would play an important 
role (e.g., Meszaros and Rees 1998; Hailey et al. 1999; Qin 
2002; Qin et al. 2004). The Doppler effect is the photons 
emitted from the region on the line of sight and those off 
the line of sight an angle of 6 are Doppler-boosted by dif- 
ferent factors and travel different distances to the observer 
Therefore, we adopt the model derived by Qin (2002) and 
Qin et al. (2004) to explore the spectral lags. But our at- 
tention is concentrated on the case of ~0. The count rate 
within energy channel [i^i, v^] is rewritten as formula (1). 



C{t) = [2nRl I{re){l + ^Tefil - r + Te)dre 

(1) 

where the tq is dimensionless relative local time defined by 
Tg = c(te — tc)/Rc, te is the emission time in the ob- 
server frame, called local time, of photons emitted from 
the concerned differential surface dsg of the fireball (6 is 
the angle to the line of sight), tc is the initial local time 



which could be assigned to any values of tg, and Rc is 
the radius of the fireball measured at tg = tc- Variable 
r is a dimensionless relative observation time defined by 
T = [c{t ~ tc) — D + Rc]/Rc, where D is the distance 
of the fireball to the observer, and t is the observation time 
measured by the distant observer. 

In formula (1), /(rg) represents the development of the 
intensity magnitude of radiation in the observer frame, called 
as a local pulse function, and ga.^{i'Q.g) describes the rest- 
frame radiation mechanisms. 

For the sake of simplicity we first adopt Gaussian pulse 
as rest-frame local pulse. As for the rest-frame radiation 
spectrum we employ the most frequently used Band func- 
tion (Band et al. 1993) since it is rather successfully em- 
ployed to fit the spectra of the GRBs. The two forms are 
rewritten as formulas (2) and (3). 
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where /o, a and Tg^min are constants, ao and /3q are the low- 
and high-energy indices in the rest frame, respectively, and 
vq^p is the rest-frame peak frequency. 

Light curves determined by formula (1) are dependent 
on the integral limits Te,min and re, max, which are deter- 
mined by the concerned area of the fireball surface, together 
with the emission ranges of the radiated frequency and the 
local time. The integral limits are only determined by 



Te,min = max{Te_r 
Tg.max = min{ri9,n: 
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where re min and rg max are the lower and upper limit of 
Tg confining I{Tg), and 6'niin and 0max are confined by the 
concerned area of the fireball surface. The radiations are ob- 
servable within the range of 



(1 - COS^min) + P COs6',„in)Te,min < 

(1 - COS ) + (1 - ^COS0max)T6( 



< T < 



.max 
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Due to the constraint to the lower limit of Tg, which is 
Tg^min > — (see, Qin et al. 2004), we assign r^.o = 
lOfT + Tg,,nin SO that the interval between Tg^min and Tg^o 
would be large enough to make the rising phase of the local 
pulse close to that of the Gaussian pulse. From formula (2) 
one can obtain 
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which leads to 
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where /S.Tg ^pw h m is the FW HM (full width at half max- 
imum) of the Gaussian pulse. From the relation between Tg 
and te, one gets 

^tg^pWHM = {Rc/c)ATg^pwHM- (8) 

In the following sections, we assign {2ttR^Iq) / {hcD^)=l, 
Rc = 3x 10^5 cm, Tg^^,„= . 

Because peak times of different light curves associated 
with different frequency intervals [(^1,(^2] or different en- 
ergy bands [Ei , E2] are different, following Shen et al. (2005) 
and Lu et al. (2006), we define the spectral lags as the dif- 
ferences between the pulse peak times in two different en- 
ergy channels, a lower energy channel [£'1, E2] and a higher 
energy channel [E^ , E4] . In the following analysis we only 
consider the energy channel pair: BATSE channels 1 (25-50 
keV) and 3 (100-300 keV) which we denote the lag as lagis 
since it was investigated by most of authors (e.g., Norris et 
al. 2000; Chen et al. 2005). 

3 Spectral lags caused by the spectral 
evolution alone 

Since Kocevski & Liang (2003) pointed out that the ob- 
served lag is the result of spectral evolution, let us first in- 
vestigate the spectral lags due to the spectral evolution alone. 
That is we consider the radiations emitted from a very small 
area with 0„iin = and 6max = 10~^ x l/F where the pho- 
tons emitted travel almost same distances to the observer In 
this case we think that the curvature effect has no role on 
the producing of spectral lag. 

Observations show that the spectral parameters of the 
GRBs vary with time. In addition, the parameters of the 
Band function low- and high-energy indices a, (3 and the 
peak energy of GRB's i^F^ spectrum Ep in the observed 

frame are mainly distributed within —1.5 0.5, —3.2 

2.1 and 100 — 700 keV and their typical values are 

a = -1 and /3 = -2.25, Ep = 250 keV, respectively 
(Preece et al. 1998, 2000; Kaneko et al. 2006). Qin (2002) 
showed the Doppler effect can not affect the intrinsic spec- 
tral shapes. When we take the Doppler effect of fireballs 
into account, the observed peak energy would be related to 
the peak energy of the rest-frame Band function spectrum 
(ao = a, /3o = ^) by Ep ~ l.eYTEo^p (see. Table 4 in 
Qin 2002). In the following we first consider the common 
evolutionary model of the rest-frame spectral parameters to 
investigate the spectral lags. Then the other possible cases 
are also taken into account. 

3.1 Spectral lags of the rest-frame spectral 
parameters remain constant 

We first explore if the spectral lags would occur when the 
rest-frame spectral parameters remain unchanged to check 
the effect of spectral evolution on the spectral lag. The typi- 
cal values of the low- and high-energy indices = — 1 and 
l3o = -2.25andi;o,p = 250/1.67rkeV are taken. Here we 



assign the a = 1.0 and take T = 100, 200, 300, 400, 500, 
600, 700, 800, 900, 1000, respectively since the Lorentz fac- 
tor of fireball are generally larger than 100. We find that all 
the spectral lags ^0913 in this case are 0. This suggests that 
there has no spectral lags when the rest-frame spectral pa- 
rameters are constant and the curvature effect does not play 
a role on the spectral lag. Fig. 1 (panel (a)) gives a exam- 
ple plot for this case. The lines of the peak times of the two 
pulses are superposed each other. 

3.2 Spectral lags of the rest-frame spectral 
parameters varying from hard to soft 

Many investigations of the spectra of GRB show that the 
spectral softening is an universal phenomenon. In addition, 
Koceviski and Liang (2003) pointed out as the Ep of lyFi, 
spectra decays through the four BATSE channels and pro- 
duce what we measure as lag. Schaefer (2004) found that 
only the evolution of Ep can produce observed lag using 
the general Liang-Kargatis relation (which describes how 
the peak photon energy in the spectrum changes with time). 
Therefore, we first investigate the case of Ep in the rest- 
frame varying from hard to soft. Following Qin et al. (2005), 
we assume a simple evolution of the peak energy -Eo.p, which 
follows logEo^p = 0.1 - k{Tg - Tgs)/{Te,2 - Te^i) (keV) 
for Tg,i <Tg< re, 2- For Tg < Tgj, log Eo,p = 0.1 (keV), 
while for Tg > Tg^2, logi?o.p — 0-1 ~ k (keV). 

Here we also take F = 100, 200, 300, 400, 500, 600, 700, 
800, 900, 1000, respectively and = -1, /^o = -2.25, 
<T ~ 1.0 are assigned. For each F we assume k =0.1, 0.2, 
0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, 1.0, 1.1, 1.2, 1.3 and 1.4 (they 
correspond to different rates of decreasing) to investigate 
the spectral lag's dependence. 

The local Gaussian pulse (see equation (2)) is employed 
to study this issue. We adopt Tgj = 8cr + Tg,min, Tg,2 = 
I2a + Tg,„iin and Tg^o = lOcr + re,™j„. 

In this way, we can compute the spectral lags. It is found 
that all of the lags are positive. A example light curve is 
showed in Fig. 1 (panel (b)), which illustrates the connec- 
tion between the spectral evolution and the spectral lag. This 
indicates the maximum of the photon count rates at the higher 
energy arrives earlier than at the lower energy in the case 
of hard-to-soft. Therefore, we can deduce that only hard- 
to-soft spectral evolution can produce the positive spectral 
lags. 

3.3 Spectral lags of the rest-frame radiation form 
varying from soft to hard with time 

Band (1997) found that a small fraction of GRB spectra 
10%) vary from soft to hard with time. Kargatis et al. (1994) 
also found this case. Therefore, we also investigate the case 
of the rest-frame radiation form varying from soft to hard. 
We also assume a simple case that the evolution of peak en- 
ergy £'o,p follows: log En^p = -0-8 + k{Tg-Tg^i) /{Tg ^2- 
Tg^i) (keV) for Tg^i < Tg < Tg-2. For Tg < Tg^i, and 
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Fig. 1 Example light curves of the Band spectral parameters Eq p remain unchanged (a), the Eq p varying from hard 
to soft (b) and soft to hard (c), where the solid and dashed curves represent the light curves of BATSE channels 1 and 3, 
respectively. 



log Eo^p = -0.8 (keV), while for Tg > Tg^2, log-Eo,p — 
—0.8 + k (keV). The F, k, and /3o are taken the same 
values as the case of hard-to-soft. We also adopt Tg i — 

So- + Tg^,nin, Tg ^2 = 12cr + Tg ^min andrg^o = lOo + Tg^jnin- 

According to our computation the soft-to-hard spectral 
evolution can result in the negative spectral lags alone. Fig. 
1 (panel (c)) also demonstrates two example light curves in- 
dicating the connection between the spectral evolution and 
the spectral lag in the case of the rest-frame radiation form 
varying from soft to hard with time. 



tween kl and k2 coming to at least 0.5. That is if the increas- 
ing rate of soft to hard is larger than the decreasing rate of 
hard to soft by at least 0.5 the negative lags would occur. As 
the increasing of F the difference needed would decrease. 
The above assumption is reasonable for the tracking pulse 
because it is consistent with the result of the hardness evo- 
lutionary characteristics of the tracking pulses investigated 
by Peng et al. (2009a) that the duration of the rise phase are 
much shorter than that of the decay phase. Fig. 2 illustrates 
four example light curves. 



3.4 Spectral lags of the rest-frame radiation form 
varying from soft to hard then to soft 

As several authors pointed out that the evolution of the low- 
energy index a and Ep exhibit the "tracking" (the observed 
spectral parameters follow the same pattern as the flux or 
count rate time profile, i.e. the evolution of spectral param- 
eters exhibit soft-to-hard-to-soft) behaviors (e.g., Ford et al. 
1995; Crider et al. 1997; and Kaneko et al. 2006; Peng et 
al. 2009a, 2009b). So we investigate the case of Eq p vary- 
ing from soft to hard then to soft. In order to ensure that 
the three parameters are in the range of observed values 
we adopted the following evolutionary form: logi?o p = 

-0.8 + kl{Tg - Tg_i)/{Tgfl - Tgs) (keV) for Te,i <Tg< 
T0,o, logi^o.p = -0.8 + fcl - k2{Tg ~ Tgfi)/{Tg^2 - Tg^) 

{keV) for t^.q < Tg < Te,2, log-Bo,p = -0.8 + fcl - fc2 
(fceV) for Tg > Tg^2, while for Tg < Tg^i, logEo^p = —0.8 
(keV). where the kl and k2 denote the increasing rate of 
soft-to-hard phase and the decreasing rate of hard-to-soft 
phase, respectively. We take F = 100, 200, 300, 400, 500, 
600, 700, 800, 900, 1000, respectively and fcl = 0.3, 0.4, 
0.5, 0.6, 0.7, 0.8, 0.9, 1.0, and 1.1, fc2 = 0.1, 0.2, 0.3, 0.4, 
0.5, 0.6, 0.7, and 0.8, respectively. We also adopt Tg^i ^ 

8ct + Tg^rnin, Tg ^2 = 12cr + Tg^.„iin SndTg^ = 10(7 + Tg^rnin- 

We use this form to investigate the spectral lags and 
find that the negative lag disappears. We explore i?o.p and 
ao varying from soft to hard then to soft. The same form 
as iJo.p is also adopted for ao. We find there are negative 
spectral lags in this case but it does not occur in the small 
F =100. When the F is larger than 100, the negative lags 
will come into being only in the case of the differences be- 



3.5 Spectral lags of a suddenly shining and gradually 
dimming intrinsic emission profile 

If the different intrinsic emission profiles are different in 
producing the spectral lags. Similar to Qin et al. (2004) and 
Shen et al. (2005) we also consider the one-sided exponen- 
tial decay and power-law decay profile, which are listed as 
follows, respectively: 

I{Tg)^hGM-C-^—^^)] iTe,n.in <Tg), (9) 



I{Tg) = Io{l - " ' r {roMn < Tg), 

'^9. max. '^0,min. 

(10) 

We assume the rg^min = and /i = 2 in the follow- 
ing analysis. The hard-to-soft and soft-to-hard spectral evo- 
lution are investigated to check if the two sorts of local 
pulses also result in the same results as those of the Gaus- 
sian pulses. The same evolutional patterns adopted above 
are employed to explore this issue. We find the results are in 
agreement with that of Gaussian local pulses. The example 
plots in Figs. 3 and 4 also illustrate the connection between 
the spectral evolution and spectral lag for the exponential 
decay and power-law decay, respectively. This suggests the 
intrinsic emission profile has no effect on the connection 
between the spectral evolution and the spectral lags. Like- 
wise, the hard-to-soft evolution produces positive lags and 
the soft-to-hard evolution leads to negative lags if we adopt 
other rest-frame local pulses. 
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Fig. 2 Example light curves of the negative (left panel) and positive (right panel) spectral lags for the case of soft-to- 
hard-to-soft, where the symbols are the same as those adopted in Fig. 1. 
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Fig. 3 Example light curve of Eo.p varies from hard to soft (left panel) and soft to hard (right panel) for the case of the 
exponential decay local pulse, where the symbols are the same as those adopted in Fig. 1. 



4 Spectral lags caused by the spectral 
evolution and the curvature effect 

As previous section pointed out that the curvature effect 
must be play a role on the producing of temporal and spec- 
tral profile we also take the effect into account to study the 
spectral lag. The evolutionary form, hard-to-soft and soft- 
to-hard, are combined with curvature effect to explore the 
influence of the combined effects on the spectral lags. That 
is we consider the radiation emitted from a big area with 
Gmin = and dmax = 'I'/S whcrc the photons emitted 
travel different distances to the observer. The photons on 
the line of sight arrive at us earlier than that of off sight- 
line. In this way the lags due to curvature effect is always 
positive and the hard-to-soft evolution also make the lag is 
positive. Therefore we think the spectral lags caused by the 
combined effect must be increase. 

We first consider the case of Eo^p varying from hard 
to soft combined with the curvature effect. In addition, we 
also assume the evolution of E'o.p follows the same means 
as the above section. To compare the lags produced by the 
spectral evolution alone with that by the spectral evolution 
and curvature effect together we plot two sorts of lags to- 



gether in Fig. 5 (panel (a)) by considering a set of parame- 
ter values. Note that the histogram plots may mean nothing 
because when one consider another set of parameter val- 
ues they get a different result. It is shown that (i) the aver- 
age value of spectral lags caused by spectral evolution and 
combined effect are 0. 1086 s and 0. 1692 s, respectively and 
(ii) the corresponding medians are 0.007 s and 0.016 s, re- 
spectively. These show the lags caused by the combined 
effect increase indeed. Whether the curvature effect also 
make the lags increase in the case of soft-to-hard or not? 
We also check the spectral lags resulting from the curva- 
ture effect and evolution of E'q p following soft-to-hard. The 
lags produced by spectral evolution alone and the combined 
effect are also compared. The mean values are -0.1200 s 
and -0.0114 s for spectral evolution and combined effect, 
respectively. Whereas the corresponding median values are 
-0.0102 s and -0.0005 s, respectively. Fig. 5 (panel (b)) also 
shows that the lags indeed increase in the case of soft-to- 
hard. In addition, we also find due to the curvature effect 
some negative lags change into positive ones. In view of the 
analysis of the above two cases we can conclude that the 
curvature effect indeed make the lags increase significantly. 
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Fig. 4 Example light curve of Eo,p varies from hard to soft (left panel) and soft to hard (right panel) for the case of power 
law decay local pulses, where the symbols are the same as those adopted in Fig. 1. 
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Fig. 5 The comparisons of lag distribution between the case of hard-to-soft (panel a) and soft-to-hard (panel b) evolution 
alone (solid lines) and the case of spectral evolution and curvature effect simultaneously (the dashed lines). 



5 Observed spectral lags of the hard-to-soft 
pulses 

The above analysis shows that the spectral lags resulting 
from the hard-to-soft spectral evolution are positive and the 
negative lags may come from the soft-to-hard or soft-to- 
hard-to-soft spectral evolution. We wonder if the observa- 
tions are consistent with the predictions. The pulses exhibit 
hard-to-soft, soft-to-hard and soft-to-hard-to-soft spectral evo- 
lution should be checked. However, the spectral parame- 
ter evolution of observed pulses generally exhibit hard-to- 
soft and "tracking" (i.e. soft-to-hard-to-soft). Both behav- 
iors are defined by Ep and/or low-energy index evolution, 
which are found by several authors (Kaneko et al. 2006; 
Ford et al. 1995; Crider et al. 1997). With ordinary Poisson 
noise in GRB light curves, there will be significant scatter in 
the measured lags. The percentage of lags increase substan- 
tially towards small values, with many (the high luminos- 
ity events) having near-zero lags. In such a case, we expect 
there to be many "measured" negative lags, even if all lags 
are positive-definite. That is, a convolution of a positive- 
definite lag distribution with the known measurement errors 



apparently reproduces the observed lag distribution includ- 
ing the measured negative lags. With this, a good case can 
be made that negative lags do not exist. The existence of 
multiple overlapping pulses could easily produce negative 
lags, even though the lags for each pulse are positive. All it 
takes is two nearly overlapping pulses with the second one 
being harder than the first. In all, we are not convinced in 
the significant existence of negative lags. 

Therefore, we select a sample only containing the FRED 
(fast rise and exponential decay) pulses provided by Ko- 
cevski et al. (2003) to check the theoretical predication since 
these bursts exhibit clean, single-peaked or well-separated 
in multi-peaked events. We can obtain directly the spec- 
tral parameters of bright bursts from Kaneko et al. (2006). 
The weaker bursts are also analyzed using RMFIT software 
and fitted with Band model. Similar to Peng et al. (2009a, 
2009b, 2010) and Ryde & Svensson (2002) a signal-to-noise 
ratio (S/N) of the observations of at least > 30 to get higher 
time resolution are adopted. For the hard-to-soft pulses only 
that containing at least 1 data points of Ep before the peak 
time of pulse and 4 in all are included in the our analysis. 
Finally, we select 15 hard-to-soft pulses as our sample to 
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Fig. 6 Spectral evolution of Ep and temporal behaviors of GRB pulses of BATSE channels 3 and 1 in our sample for the 
case of hard-to-soft, where the black lines represent the hght curves of channel 3. The positive lags are evident. 
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test our predictions. Illustrated in Fig. 6 are the spectral and 
temporal behaviors for our sample. 

We first employ the widely used cross-correlation func- 
tion (CCF) to compute the spectral lags between the pulses 
of the same burst seen in the BATSE channels 1 and 3. Sim- 
ilar to Norris et al. (2000) we adopt a cubic polynomial to 
fit the peak of the resulting discrete CCF function and take 
the peak of the cubic polynomial form as the lag. The error 
in the average CCF lag, lagis, is found by simulations with 
Gaussian noise added to each of the two channels. These 
spectral lags resulting from the evolution of hard-to-soft are 
listed in Table 1 . We find no negative lags in the hard-to-soft 
pulses. 

Recendy Hakkila et al. (2008) point out that the CCF 
and pulse peak lags occasionally disagree. In order to fur- 
ther test it we also find pulse peak lags which are defined 
as the differences between the pulse peak in BATSE chan- 
nels 1 and 3. In order to find the peak time of pulse we must 
select a pulse model to fit these pulses. Following Peng et 
al. (2006) we only select the function presented in equa- 
tion (22) of Kocevski et al. (2003). We fit those selected 
background-subtracted pulses using the KRL function. Our 
fittings are constrained on the channels 1 and 3 to find their 
difference of peak times. The best fitting parameters are de- 
pended on which xt of the fitting model is smaller The 
pulse peak lags along with the fitting xf, in channels 1 and 3 
are also listed in Table 1 . From Table 1 we find all the lags 
of hard-to-soft pulse are positive. We also check the CCF 
lags and peak lags by smoothing the light curves with the 
DB3 wavelet (Qin et al. 2004) and obtain the same results 
(see. Table 1). 

6 Conclusions and discussions 

Spectral evolution is an established characteristic in GRBs. 
Moreover the trend for their high-energy photons to arrive 
before the lower-energy ones ("hard-to-soft" evolution) is 
universal phenomenon, which leads to positive lags. Ko- 
cevski & Liang (2003) showed the robust connection be- 
tween the spectral lag measured in GRBs and the hard-to- 
soft evolution of the burst spectra by employing a sample of 
clean single-peaked bursts with measured lag. In this paper, 
we have investigated the connection between the evolution 
of the rest-frame spectral parameters and the spectral lag 
using a theoretical model derived by Qin (2002) and Qin 
et al. (2004). We first assume an intrinsic Band spectrum 
and a Gaussian emission profile to study the issue. Our the- 
oretical analysis shows that the connection exists not only 
between the hard-to-soft evolution of the burst spectra and 
the spectral lag but between the soft-to-hard evolution and 
the spectral lag as well as between the soft-to-hard-to-soft 
and the spectral lag. In addition, the hard-to-soft spectral pa- 
rameters evolution can produce the positive lags, while the 
negative lags may be resulting from the evolution of soft- 
to-hard as well as soft-to-hard-to-soft. We also take account 
of the other local pulse forms studied by Qin et al. (2004) 
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Fig. 7 Peak lag (s) versus CCF lag (s). A good linear cor- 
respondence exists, as expected. 

and find that the spectral lags resulting from these different 
local pulses are consistent with those of the Gaussian pulse. 

Let us check whether there are any differences in the 
measurements between peak lags, Lagi^, adopted in the 
above analysis and the widely used cross correlation func- 
tion lags, CCFi3 (e.g., Norris et al. 2000). Similar to Band 
(1997) and Norris et al. (2000) we fitted a cubic polyno- 
mial to the peak of the resulting discrete CCF function and 
find the difference of the corresponding peak time. In Fig. 
7 the peak lag, Lagis, is plotted as a function of the CCF 
lag, CCF13. A linear fit is given which shows that there is a 
good correspondence between the two methods. The best fit 
is Lagi3 = (-0.03 ± 0.03) + (1.28 ± 0.23) x CCF13. The 
values of intercept and slope are consistent with the results 
of Ryde et al. (2005). Therefore, the peak lags we used in 
this study would not affect our analysis results. 

We select a sample consisting of 15 pulses whose peak 
energy Ep follow hard-to-soft evolution to test our theoreti- 
cal predictions. Our analysis shows that all of the CCF lags 
and analytic lags of hard-to-soft pulses are positive within 
the range of uncertainty. In order to decrease the effect of 
noise on our results we smoothed the data with the DB3 
wavelet (the second-class decomposition) using the MAT- 
LAB software (Qin et al. 2004). Then we re-calculate the 
two types of lags (see. Table 1) and find the lags are also 
positive. This indicates hard-to-soft spectral evolution may 
indeed produce positive lags, which is in good agreement 
with our theoretical prediction. 

Currently we can not ascertain if the negative lag pre- 
dicted by our adopted model are consistent with observa- 
tion because we are not convinced in the significant exis- 
tence of negative lags in our observation. The main reasons 
mentioned in the previous section are the near-zero lags of 
most pulse as well as the existence of multiple overlapping 
pulses. Even then we can make a theoretical predication 
about the producing of the negative lags, i.e. the soft-to-hard 
and the soft-to-hard-to-soft spectral evolution may give rise 
to the negative lags, which reserves the further investiga- 
tions in the future observations. 
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Table 1 A list of various observed lags of GRB pulses whose spectral evolution following from hard to soft 



trigger 


St (s) 


et (s) 


iagrs,p (s) 


dlagis^p (s) 


iagia,c (s) 


dlagis,c (s) 




xi,3 


1406 


-2.0 


15.0 


1.864 ±0.169 


1.857 ±0.169 


1.152 ±0.019 


0.487 ± 0.509 


1.05 


1.07 


2138 


-2.0 


15.0 


0.988 ± 0.302 


1.001 ±0.316 


1.024 ± 0.057 


0.604 ±0.178 


1.10 


0.99 


2387 


-2.0 


20.0 


2.188 ±0.204 


2.169 ±0.207 


1.856 ±0.019 


1.574 ±0.177 


1.12 


1.15 


2662 


-2.0 


15.0 


1.331 ±0.285 


1.355 ±0.272 


0.768 ±0.159 


0.712 ±0.139 


1.01 


1.17 


2665 


-2.0 


15.0 


1.552 ±0.196 


1.555 ±0.196 


1.507 ±0.503 


1.417 ±0.148 


0.89 


0.91 


3256 


-2.0 


15.0 


2.120 ±0.428 


2.828 ± 0.390 


1.532 ±0.806 


1.325 ±0.147 


1.04 


1.07 


3257 


-1.0 


15.0 


1.667 ±0.517 


1.606 ±0.509 


1.283 ±0.075 


1.134 ±0.094 


0.97 


0.92 


3875 


-1.0 


5.0 


0.262 ± 0.029 


0.257 ±0.031 


0.162 ±0.026 


0.169 ±0.011 


1.40 


1.18 


4368 


-2.0 


10.0 


0.346 ± 0.047 


0.346 ± 0.048 


0.265 ± 0.027 


0.265 ± 0.013 


1.03 


0.82 


5621 


3.2 


4.5 


0.032 ± 0.022 


0.076 ±0.011 


0.028 ± 0.002 


0.030 ± 0.001 


1.24 


1.71 


6397 


-2.0 


20.0 


1.158 ±0.082 


1.160 ±0.082 


0.942 ±0.031 


0.934 ±0.171 


1.17 


1.70 


7548 


2.0 


12.0 


0.814 ± 0.075 


0.816 ±0.075 


0.596 ± 0.232 


0.628 ± 0.063 


0.88 


1.10 


7588 


-2.0 


15.0 


2.230 ±0.120 


2.220 ±0.120 


1.125 ±0.292 


1.427 ±0.390 


0.87 


0.93 


7638 


-2.0 


15.0 


1.124 ±0.114 


1.420 ±0.134 


0.968 ± 0.012 


0.931 ±0.162 


1.15 


1.04 


7648 


-2.0 


15.0 


2.490 ± 0.243 


2.501 ± 0.248 


1.574 ± 1.202 


1.528 ±0.386 


1.01 


1.03 



Note: St and et denotes the start and end time since bursts trigger of selected pulses, respectively, lagis^p and lagis^c represents peak 
time lag and CCF lag, while dlagri,p and dlagn.c are the corresponding lags but the light curve having been smoothed with the DB3 

wavelet, and x^,3 is the fitting xt of channels 1 and 3, respectively. 



When taking the spectral evolution and curvature ef- 
fect into account simultaneously we find the lags are larger 
than those associated with the case that we only consider 
the spectral evolution (see. Fig. 5). The two factors must 
play the roles simultaneously on the producing of the spec- 
tral lags since the curvature effect must be play a important 
role on the producing of spectral and temporal profile of 
the GRB. In other words, the observed lags should suffer 
from this two factors simultaneously. Therefore we argue 
that the number of the observed negative lag should be much 
smaller than positive lags. Recently, Chen et al. (2005) used 
the data acquired in the time-to-spill (TTS) mode for long 
GRBs by BATSE and found that positive lags is the norm 
in this selected sample set of long GRBs. While relatively 
few in number, some pulses of several long GRBs do show 
negative lags, which is also consistent with our theoretical 
predication. 

Both the theory and observation suggest the connection 
between the spectral evolution and the spectral lags indeed 
exists. Therefore, we must make clear the mechanism of in- 
trinsic spectral evolution in order to reveal the underlying 
physics of the spectral lag. To explore this, we must first 
understand the mechanism that produces breaks in the GRB 
spectra and what causes its evolution but it is a difficult is- 
sue. As for the physical model of hard-to-soft spectral evo- 
lution several authors have proposed the possible origin. For 
example, Tavani ( 1 996) showed the behavior of hard-to-soft 
spectral evolution are caused by the variation of the average 
Lorentz factor of pre-accelerated particles and the strength 
of the local magnetic field at the GRB site as the synchrotron 
emission evolves within the burst. While Liang (1997) pro- 
posed a physical model of hard-to-soft spectral evolution in 
which impulsively accelerated non-thermal leptons cool by 
saturated Compton upscattering of soft photons. Kocevski 
& Liang (2003) pointed out it is a harder question due to un- 
certainties involved with the microphysics of GRBs. They 



thought the interpretation of peak energy Ep (and hence its 
evolution) depends on the radiation mechanism that is used 
to explain the GRB spectra. Moreover they also attempted to 
examine several mechanisms that could produce the decay 
of the GRB spectra and discuss how this evolution could be 
connected to the bursts luminosity for each model. Whereas 
for the soft-to-hard-to-soft spectral evolution few attempts 
have been made to explain it. Peng et al. (2009a) argued 
that the phenomenon may be caused by both kinematic and 
dynamic process. However they did not provide further ev- 
idences to distinguish which one is dominant. The issue of 
spectral evolution of hard-to-soft and soft-to-hard-to-soft is 
still unclear, which deserves careful studies and further in- 
vestigations. 
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